Are Light Sterile Neutrinos Preferred or Disfavored by Cosmology? 
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We find that the viability of a cosmological model that incorporates 2 sterile neutrinos with 
masses around 1 eV each, as favored by global neutrino oscillation analyses including short baseline 
results, is significantly dependent on the choice of datasets included in the analysis and the ability 
to control the systematic uncertainties associated with these datasets. Our analysis includes a 
variety of cosmological probes including the cosmic microwave background (WMAP7+SPT), Hubble 
constant (HST), galaxy power spectrum (SDSS-DR7), and supernova distances (SDSS and Union2 
compilations). In the joint observational analysis, our sterile neutrino model is equally favored as a 
ACDM model when using the MLCS light curve fitter for the supernova measurements, and strongly 
disfavored by the data at Axoft ~ 18 when using the SALT2 fitter. When excluding the supernova 
measurements, the sterile neutrino model is disfavored by the other datasets at Ax^a ~ 12, and 
at best becomes mildly disfavored at A-XeS ~ 3 when allowing for curvature, evolving dark energy, 
additional relativistic species, running of the spectral index, and freedom in the primordial helium 
abundance. No single additional parameter accounts for most of this effect. Therefore, if laboratory 
experiments continue to favor a scenario with roughly eV mass sterile neutrinos, and if this becomes 
decisively disfavored by cosmology, then a more exotic cosmological model than explored here may 
become necessary. 



I. INTRODUCTION 

The standard models of particle physics and cos- 
mology do not yet fully describe the neutrino sector, 
with open questions related to the mass-generation 
mechanism of the neutrinos, any sterile neutrino part- 
ners of the active neutrinos, and their potential rela- 
tion to the number of relativistic degrees of freedom 
inferred from cosmology. In recent years, there has 
been some experimental evidence pointing towards 
the existence of additional light (effectively masslcss) 
degrees of freedom. In particular, a combined anal- 
ysis of cosmic microwave background (CMB) data 
from WMAP7, baryon acoustic oscillation (BAO) dis- 
tances from SDSS+2dF, and Hubble constant from 
HST yields a weak preference for additional light de- 
grees of freedom (N cS = 4.34 ±0.87) Q. When more- 
over including small-scale CMB data from ACT or 
SPT, this preference mildly increases to the 2a level 
(N eS = 4.56 ± 0.75 with addition of ACT @ and 
N cS = 3.86 ± 0.42 with addition of SPtJI). This 
possibility has sparked further work |J^2CJ. These 
constraints on N c ff explicitly assume that the addi- 
tional particles are masslcss. 

However, in light of new predictions for the anti- 
neutrino flux from nuclear reactors, global short- 
baseline neutrino oscillation data now favor the ex- 
istence of two sterile neutrinos with best-fit masses of 
TO4 = 0.68 eV and m; = 0.94 eV, assuming mass- 
less active neutrinos [5l| (also see [22|-[25j]). Instead 
of analyzing the data with the aim of estimating an 
upper bound to the mass of an additional thcrmalizcd 
neutrino species [H, [TH, [l7[ , we take the existence of 
two sterile neutrinos with m& and m$ as a prior as- 
sumption consistent with the short-baseline data. It 
is our aim to determine how a model with these two 
additional neutrino species fares compared to the case 



Parameter 


Symbol 


Prior 


Baryon density 


n b h 2 


0.005 -» 0.1 


Cold dark matter density 


n c h 2 


0.01 -> 0.99 


Angular size of sound horizon 


e s 


0.5 -» 10 


Optical depth to reionization 


r 


0.01 0.8 


Scalar spectral index 


n 3 


0.5 -> 1.5 


Amplitude of scalar spectrum 


In (lO 10 ^) 


2.7 -> 4 


Effective number of neutrinos 


iV ea 


3.046 10 


- with sterile neutrinos 


N cH 


5.046 -»■ 10 


Sum of neutrino masses 


[eV] 





- with sterile neutrinos 


E™* [eV] 


1.62 


Constant dark energy EOS 


w 


-3^0 


Running of the spectral index 


dn s 
d In k 


-0.4 -> 0.4 


Curvature of the universe 




-0.4 -»■ 0.4 


Primordial helium abundance 


Y P 


1 



TABbE I. We impose uniform priors on the above cosmo- 
logical parameters. In addition, we always consider the 
Poisson point source power D^oo, the clustered power 
D3Q00, and the SZ power -Dfoon as nuisance parameters 
constrained by the CMB data [g]. Moreover, we always 
derive as, the amplitude of linear matter fluctuations on 
scales of 8 Mpc/h at z = 0. We only vary a redshift- 
independent dark energy equation of state (EOS). In this 
table, the first 6 parameters are defined as "vanilla" pa- 
rameters. 



without them, when including all available and rele- 
vant cosmological data. 

We examine the impact of the two sterile neutri- 
nos on other cosmological parameters in the vanilla 
ACDM model, such as the matter density, amplitude 
of linear matter fluctuations on 8 Mpc/h scales, and 
spectral index. We also explore the impact of ex- 
tensions of a cosmological model with sterile neutri- 
nos, including nonzero curvature, evolving dark en- 
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FIG. 1. Left: CMB temperature power spectrum measurements with WMAP7 (orange) and SPT (blue). The ACDM 
model without sterile neutrinos is shown with the solid (black) line, and the ACDM model with 2 sterile neutrinos is 
shown in dashed (red). Right: Assuming the ACDM model is centered on the DR7 data, with error bars given by the 
shaded band (in blue), we show the absolute difference with our sterile neutrino model in solid (red). 



ergy, running of the spectral index, and primordial 
helium abundance. Throughout this paper, we will 
assume that the two sterile neutrinos are thermally 
populated. If this is not the case, then the differences 
between a model with two sterile neutrinos and one 
without them will be smaller (cf. Ref. (26j). 

The cosmological influence of sterile neutrinos in- 
cludes an increase in the effective number of neutri- 
nos to A^off = 5.046 and the sum of neutrino masses 
to yimy = 1.62 eV, since the sterile neutrinos mix- 
ing angles would require their thermalization in the 
early universe [13, [H| (see, however, Ref. [26| for 
cases where only partial thermalization may occur). 
As discussed in Ref. [3| , the effective number of neu- 
trinos is mainly correlated with the matter density 
and spectral index in a vanilla ACDM model. In 
extended cosmological models, correlations also exist 
with the helium abundance, dark energy equation of 
state, and running of the spectral index. Meanwhile, 
the sum of neutrino masses is mainly correlated with 
the matter density and Hubble constant in a vanilla 
ACDM model, along with the dark energy equation 
of state and curvature density in extended parameter 
spaces 

The radiation content of the universe can be 
constrained from big bang nucleosynthesis (BBN) 
through its effect on the expansion rate (29l - [3l1 | . Given 
the standard BBN consistency relation between the 
set of parameters {Y p , N c s , VL^h 2 } [30| . the inclusion 
of 2 additional neutrinos boosts the primordial he- 
lium abundance by AY p = 0.024 when the baryon 
density is kept fixed. Thus, Y p w 0.27 in standard 
cosmological analyses when enforcing this consistency 
relation. Primordial helium abundance estimations 
from observations of metal poor cxtragalactic H II 
regions suffer from significant systematic uncertain- 
ties (e.g. see [32M371 ]). An extensive analysis that at- 
tempts to account for these systematic uncertainties 
gives Y p = 0.2534 ± 0.0083 [37|, which is consistent 



with the cosmological estimate at 95% CL (assuming 
5 light neutrinos). This agreement could be tightened 
by lowering Y p from cosmology, achieved via mecha- 
nisms such as incomplete thermalization, presence of 
a non-zero chemical potential, or post-BBN produc- 
tion of the sterile neutrinos from the decay of a heavy 
particle species (e.g. see [HI). 

We describe our analysis method in Section 2. In 
Section 3, we provide constraints on a ACDM model 
with three massless active neutrinos and two massive 
sterile neutrinos, then follow up by allowing for evolv- 
ing dark energy, universal curvature, running of the 
spectral index, additional rclativistic species, and free- 
dom in the primordial helium abundance (all param- 
eters defined in Table []}. Section 4 concludes with a 
discussion of our findings. 



II. METHODOLOGY 

We employed a modified version of CosmoMC [H, 
[39( ] in performing Markov Chain Monte Carlo 
(MCMC) analyses of parameter spaces with sterile 
neutrinos, using CMB data from WMAP7 Q an d 
SPT luminous red galaxy power spectrum mea- 
surements from SDSS DR7 [40|, the Hubble con- 
stant from HST Q , and SN distances from either the 
Union2 compilation [4l[ or the SDSS compilation [42[ . 
We generally impose a cutoff in the galaxy power spec- 
trum measurements at k = 0.1 /i/Mpc because of in- 
sufficient understanding of the matter power spectrum 
on nonlinear scales when including baryons, massive 
neutrinos, and dark energy (43T4521 ] . For the same rea- 
sons, we do not include the small-scale power spec- 
trum from Lyman-a forest data. 

The Union2 compilation consists of 557 SNc, which 
includes large samples from SCP, SNLS, ESSENCE, 
HST, and older data sets gl|, while the SDSS 
compilation consists of 288 SNe from SDSS, SNLS, 
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ESSENCE, HST, and a set of low-rcdshift SNe 0. 
For the Union2 compilation, we considered the SALT2 
light curve fitter |54| , while for the SDSS compilation, 
we considered both the MLCS 0E1 and SALT2 fit- 
ters. The two fitting methods make different assump- 
tions about the nature of color variations in type la 
SNe and employ different ways in determining model 
parameters (for further details, see Ref. j42[). At 
present, there seems to exist no consensus on which 
light curve fitter is the most accurate (e.g. 0, 

All parameters are defined in Table HI The power 
spectra of the CMB temperature and E-mode polar- 
ization were obtained from a modified version of the 
Boltzmann code CAMB [57l[58j. We used the Gelman 
and Rubin R statistic [59( to determine the conver- 
gence of our chains, where R is defined as the vari- 
ance of chain means divided by the mean of chain 
variances. In stopping the runs, we generally required 
the conservative limit (R — 1) < I0~ 2 , and checked 
that further exploration of the tails does not change 
our results. 

In our baseline ACDM model, we include 3 mass- 
less neutrinos. We also consider an expanded ACDM 
model that contains 2 sterile neutrinos in addition to 
the 3 active neutrinos of the baseline model. The 
sterile neutrino masses are given by the mass split- 
tings with the lightest neutrino mass: vhA = 0.68 eV 
and TO5 = 0.94 eV [2lJ. Beyond the 3 massless ac- 
tive species and 2 massive sterile species, additional 
contributions to N e g are assumed massless. 

For the primordial fraction of baryonic mass in he- 
lium, there are three reasonable priors we can explore: 
I) fixing Y p to a constant, 2) allowing Y p to vary as 
a free parameter, and 3) determining Y p as a func- 
tion of {N e g, fife/i 2 } in a manner consistent with BBN 
(e.g. see Eqn 1 in Ref. Q). We show results when fix- 
ing the the primordial helium abundance to the SPT 
preferred value of Y p = 0.2478 0]. We have checked 
that our results do not significantly vary when forcing 
Y p to preserve the standard BBN consistency relation 
instead. As part of our analysis of extended parameter 
spaces, we also consider cases with the helium abun- 
dance as an unknown parameter to be determined by 
the data. 

We define the running of the spectral index 
dn s /dhik through the dimensionless power spectrum 
of primordial curvature perturbations: 

(/ \ n s —l-\-^ \n(k/ ko)dn s / d\n k 
k~o) ' (1) 

where the pivot scale ko = 0.002/Mpc. Due to the 
large correlation between n s and dn s /d\nk at this 
scale, we consistently quote our values for n s at a 
scale fco = 0.0I5/Mpc, where the tilt and running 
are less correlated, such that n s (fco = 0.015/Mpc) = 
n s (k = 0.002/Mpc) + ln(0.015/0.002)dn s /dlnfc [Hq|. 
An example of the remaining correlation between the 
spectral index and its running is shown in Ref. 
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FIG. 2. Joint two-dimensional marginalized constraints 
on o"8(f2m/0.25) ' 47 against €l m from combining the mea- 
surements of WMAP+SPT+P(fc)+HST+SNe. The pur- 
ple and pink shaded confidence regions (inner 68%, outer 
95%) are obtained using SNe from the Union2 compila- 
tion (SALT2), while the solid and dashed transparent el- 
lipses are obtained using SNe from the SDSS compilation 
(MLCS). The overlapping ellipses preferring a lower mat- 
ter density {left) are for the ACDM model without sterile 
neutrinos, while the overlapping ellipses preferring a larger 
matter density {right) axe for the ACDM model with sterile 
neutrinos. The horizontal dashed lines (in brown) denote 
the 68% confidence interval about the mean from the local 
(0.025 < z < 0.25) galaxy cluster abundance measurement 
of Vikhlinin et al. (2009) [gj. 



We define xts = -21n£ max , where £ max is the 
maximum likelihood of the data given the model. The 
ratio of maximum likelihoods given two separate mod- 
els is then C max , 2 / C max ,i = cxp(-Axo ff /2). For the 
case where Axle > 0, we interpret model 2 to be as- 
sociated with a lower probability of drawing the data 
at the maximum likelihood point than model I, by a 
factor given by exp(— Ax 2 ff /2). For reference, a value 
of A^eff = 10 corresponds to odds of I in 148, which 
we take as strong preference for model 1 as compared 
to model 2. 

We also consider the Deviance Information Crite- 
rion (DIC) HU, given by DIC = x 2 cS (d) + 2C b , where 

Cb = Xeff(^) ~ Xeff(^) * s the so-called "Bayesian com- 
plexity," such that 8 is the vector of varied parame- 
ters, the bar denotes the mean over the posterior dis- 
tribution, and hat denotes the maximum likelihood 
point [62(. The Bayesian complexity can be thought 
of as the effective number of unconstrained parame- 
ters, such that it penalizes more complex models with 
more parameters, independently of how well the mod- 
els fit the data [63| . If the Bayesian complexity of two 
models is the same, the difference in DIC between 
the models matches their difference in x\s va hies. We 
take a difference beyond 10 in DIC values between two 
models to constitute a strong preference for one model 
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TABLE II. Constraints on Cosmological Parameters using SPT+WMAP+P(fc)+7io- In some of the columns, we further 
add SNe from either the Union2 or SDSS compilations. The foreground priors on the SZ, poisson point sources, and 
clustering point sources are encapsulated in "FG." 







ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 








+2u a 


±SNe Union 2 


f 2i/ s +SNe Union 2 


±SNesDSS 


±2^ s +SNe S DSS 


Primary 


100Q b h 2 


2.242 ± 0.039 


2.296 ± 0.040 


2.241 ±0.039 


2.308 ± 0.040 


2.225 ±0.038 


2.293 ±0.038 




ioon dm ?i 2 


11.15 ±0.32 


16.49 ± 0.41 


11.18 ±0.29 


16.09 ±0.36 


11.63 ± 0.29 


16.51 ±0.36 




10 4 6» a 


104.15 ±0.15 


103.86 ± 0.15 


104.15 ±0.15 


103.90 ±0.15 


104.10 ± 0.15 


103.85 ±0.15 




T 


0.086 ±0.014 


0.089 ±0.014 


0.087 ±0.014 


0.092 ± 0.015 


0.082 ±0.013 


0.089 ± 0.014 




100n 3 


96.73 ±0.95 


98.32 ± 0.97 


96.69 ± 0.94 


98.81 ±0.95 


95.98 ± 0.94 


98.27 ±0.94 




ln(10 10 A s ) 


3.187 ±0.035 


3.214 ±0.036 


3.189 ±0.035 


3.195 ± 0.036 


3.216 ±0.035 


3.216 ±0.035 


Derived 


H 


71.4 ± 1.4 


69.6 ± 1.3 


71.2 ± 1.3 


71.0 ±1.2 


69.2 ± 1.2 


69.5 ± 1.2 




CT8 0W0.25) - 47 


0.829 ± 0.039 


0.813 ±0.039 


0.833 ± 0.035 


0.776 ± 0.034 


0.886 ±0.036 


0.816 ±0.036 


V 2 

A-eff 


CMB 


7512.4 


7517.2 


7511.7 


7517.9 


7513.2 


7516.7 




P(k) 


23.9 


28.9 


24.5 


30.2 


23.2 


28.7 




Ho 


1.5 


2.9 


1.2 


1.5 


4.6 


3.5 




SNe 






530.8 


536.0 


245.9 


237.9 




FG 


0.1 


0.6 


0.4 


0.7 


0.1 


0.7 




Total 


7537.9 


7549.5 


8068.5 


8086.2 


7787.0 


7787.4 


DIC 


Total 


7554.1 


7566.1 


8085.2 


8103.0 


7803.4 


7804.1 




Total 




11.6 




17.7 




0.4 


ADIC 


Total 




12.0 




17.8 




0.7 



Mean of the posterior distribution of cosmological parameters along with the symmetric 68% confidence interval about the 
mean. The three active neutrinos are taken to be massless in all models. We also consider adding 2 sterile neutrinos (denoted as 
"2!A," ) of masses m„ sl = 0.68 and m Us2 = 0.94, such that the sum of neutrino masses is m„ = 1.62 eV. We fix the primordial 
helium mass fraction Y p = 0.2478. The Deviance Information Criterion is defined as DIC = 2x 2 f j(#) — Xjjff (^)i w h ere 6 is the 
vector of varied parameters, the bar denotes the mean over the posterior distribution, and hat denotes the maximum likelihood 
point. For the SDSS SNe, we have used the MLCS light curve fitter. The corresponding total Ax 2 ff and ADIC values when 
using the SALT2 fitter are Ax 2 ff = 20.1 and ADIC = 19.4. For the Union2 SNe, we always use the SALT2 fitter. In the rows 
with x 2 fj = — 21n£ max values listed for individual probes, the values are computed at the maximum likelihood estimate (MLE) 
for the joint analysis including all probes. If each probe is analyzed separately, the MLE will be different and the corresponding 
Ax 2 g values will be smaller. 



as compared to the second model, with the more pre- 
ferred model being the one with the smaller DIC value. 



III. RESULTS 

We now explore the cosmological constraints on our 
sterile neutrino models, and the relative goodness of 
fit with respect to models without sterile neutrinos. In 
Sec. IIII Al we vary the parameters of a vanilla model 
defined in Table [H while we consider an extended pa- 
rameter space in Sec. IIII B[ 



A. Vanilla plus 1v s Models 

In Table UH we show the constraints on two sep- 
arate ACDM models for three distinct supernova 
cases: 1) without SNe, 2) with Union2 SNe (SALT2 fit- 
ter), and 3) with SDSS SNe (MLCS fitter). The model 
denoted "ACDM" consists of the 6 vanilla parameters 
in Tablc[jjand does not contain sterile neutrinos, while 
the model denoted "ACDM+2i/ s " consists of the same 
vanilla parameters but now contains two sterile neu- 



trinos of fixed masses vn± = 0.68 eV and m§ = 0.94 eV 
(as discussed in Scc.|TT|. We define &XcS as being the 
difference in Xeg between the sterile neutrino model 
(ACDM + 2v s ) with the null model (ACDM). 

When excluding SN data, we find that the model 
with sterile neutrinos is disfavored at Axis = 11-6, 
which implies a factor of 330 larger odds for the null 
model to draw the data than the sterile neutrino 
model assuming the maximum likelihood parameters. 
Moreover, ADIC = 12.0. Since ADIC ~ A^ ff , this 
tells us that the two models have essentially the same 
Bayesian complexity, and both statistical measures 
(DIC and Xcs) strongly favor the null model over the 
one with two massive sterile neutrinos. 

Allowing for SN data, the corresponding results 
are Ax^ = 17.7 and ADIC = 17.8 for the Union2 
compilation (SALT2 fitter), while Ax^ s = 0.4 and 
ADIC = 0.7 for the SDSS compilation (MLCS fitter). 
When using SDSS SN data with the SALT2 light curve 
fitter, the corresponding results are A%^ ff = 20.1 and 
ADIC = 19.4. In other words, the choice of light curve 
fitter has a decisive impact on the statistical viability 
of the sterile neutrino model. The two fitters are also 
associated with slight differences in the inferred mat- 
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1.02 



0.98 



0.96 



0.94 



0.2 



ACDM 




1.02 



0.98 



0.96 



0.94 




0.25 



0.3 



0.35 



0.3 



0.35, 



0.4 



0.45 



m m 
FIG. 3. Joint two-dimensional marginalized constraints on the spectral index n s and matter density Q, m (inner 68%, 
outer 95%). The green shaded ellipses are for WMAP+HST, blue shaded ellipses are for WMAP+SPT+HST, the 
solid transparent ellipses (in red) are for WMAP+P(£;)+HST, and the dashed transparent ellipses (in black) are for 
WMAP+HST+SNe, where the SNe are from the SDSS compilation (MLCS). The panel to the left assumes a ACDM 
model without sterile neutrinos, while the panel to the right includes two sterile neutrinos. 



ter density and Hubble constant, with larger values 
of the former and smaller values of the latter being 
associated with the MLCS fitter. These discrepancies 
may ultimately be traced back to differences between 
the fitters in the rest-frame U-band region [4l|, [42|, [65[ . 

We note that our results using the Union2 data 
set are in agreement with those in Re f. 16611 . which 
used SNe from the Union2 compilation [67| , extended 
the galaxy power spectrum measurements out to 
k = 0.2 /i/Mpc, and excluded small-scale CMB data. 
When extending the power spectrum measurements 
from fcmax = 0.1 h/Mpc out to /c max = 0.2 h/Mpc, 
Axcft increases by about 5 for all of the different cases 
including SN data. 

In Fig.[TJ we show the CMB temperature and galaxy 
power spectra for a ACDM model without sterile neu- 
trinos and one with 2 sterile neutrinos. While the 
influence of additional neutrinos is a systematic sup- 
pression in both spectra, the figures show that this 
level of suppression largely lies within the error bars of 
present data. In other words, the figures show that the 
sterile neutrino model provides a good fit to the data, 
albeit slightly worse than the null model. In Table [H] 
we directly show from which probes the largest differ- 
ences in XcS arise for our two models. For the case 
without SNe, ^XeS receives an equal contribution of 
about 5 from each of the CMB and galaxy probes. For 
the case with SNe from the Union2 compilation, A% of j 
receives an equal contribution of about 6 from each of 
the CMB and galaxy probes, and roughly 5 from the 
SNe measurements. Hence, no single probe manages 
to decisively disfavor the sterile neutrino model. 

For the case with SNe from the SDSS compilation, 
Axpff receives a contribution of 5.5 from P(k), 3.5 
from the CMB, but then a negative contribution of 
8.0 from the SN measurements (MLCS). Thus, the 
main difference between our Union2 and SDSS super- 
nova cases (with the different light curve fitters) is that 



the former disfavors sterile neutrinos, while the latter 
prefers sterile neutrinos. We note that the individual 
Xeff values in Table |H] are those associated with the 
maximum likelihood point of the joint analysis of all 
considered probes. When each probe is analyzed sep- 
arately, the best-fit A%;f ff values are less pessimistic. 

In Table HH we further show the constraints on a 
range of cosmological parameters. In particular, we 
find that the sterile neutrino model prefers a larger 
matter density and lower value of aa , while preserving 
the constraint on <Tg,(£l m /0.25) OA7 near the 0.8-mark, 
in agreement with the galaxy cluster abundance mea- 
surement of Vikhlinin et al. (2009) jgl]. In Fig.d we 
show error ellipses in the plane of (78(^m/0.25) ' 47 and 
fl m for the case of WMAP+SPT+P(fc)+HST+SNe. 
Remarkably, when using SNe from the SDSS compila- 
tion, a much larger matter density is allowed to con- 
stitute the energy content of our universe. 

In our MCMC analyses, we minimize the total 
XcS: which is the sum of the individual Xcfi va ^" 
ues of equally weighted probes. However, given the 
different systematics, the motivation for weighting 
the CMB and large-scale structure probes equally is 
not clear. In this regard, in Fig. [31 we show er- 
ror ellipses for n s against Q m given different sets 
of probes: f ) WMAP+HST, 2) WMAP+HST+SPT, 
3) WMAP+HST+P(fc), and 4) WMAP+HST+SNe. 
We find that these different combinations of probes 
constrain a portion of parameter space in agreement 
with each other, even for the case of sterile neutrinos. 
In other words, the preferred parameter space is not 
driven by a single probe, but consistently preferred by 
all probes. 

We have shown that the possibility of a cosmolog- 
ical model that incorporates 2 sterile neutrinos with 
roughly eV masses is significantly dependent on the 
choice of datasets included in the analysis and the 
ability to control the systematic uncertainties associ- 
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TABLE III. Constraints on Cosmological Parameters using SPT+WMAP+P(fc)+_H"o- In some of the columns, we 
further add SNe from the Union2 compilation adopting the SALT2 fitter. Note that the MLCS fitter does not disfavor 
the addition of 2 light sterile neutrinos and those results are not shown here since the main motivation is to gauge the 
effect of additional cosmological parameters. 







ACDM 


ACDM 


ACDM 


ACDM 








+2u a 


+SNc Union2 


±2^+SNc Union2 


Primary 


ioon b h 2 


2.235 ±0.073 


2.232 ±0.067 


2.235 ±0.076 


2.208 ±0.069 




100H dm /i 2 


13.7 ± 1.2 


19.8 ±1.4 


12.88 ±0.93 


18.6 ± 1.0 




10 4 6» s 


103.63 ± 0.37 


103.03 ± 0.28 


103.84 ± 0.35 


103.17 ±0.29 




r 


1.002 ±0.019 


0.090 ±0.016 


0.093 ±0.016 


0.084 ±0.014 




100n s 


97.5 ±2.7 


97.9 ±2.1 


97.3 ±2.3 


96.2 ± 1.7 




ln(10 10 ^ s ) 


3.105 ±0.068 


3.183 ±0.064 


3.156 ±0.045 


3.237 ± 0.044 


Extended 


n k 


0.023 ±0.021 


0.019 ±0.022 


0.003 ±0.010 


-0.004 ±0.012 




w 


-0.76 ±0.20 


-0.80 ±0.32 


-0.999 ± 0.099 


-1.32 ± 0.16 




N e g 


4.22 ±0.74 


6.83 ±0.97 


3.77 ±0.50 


6.02 ±0.67 




dn a /d In k 


-0.048 ± 0.036 


-0.029 ± 0.027 


-0.027 ±0.033 


-0.019 ±0.027 




Y P 


0.165 ±0.084 


0.086 ±0.061 


0.214 ±0.078 


0.113 ±0.067 


Derived 


Ho 


72.6 ±2.4 


72.8 ±2.5 


73.6 ±2.2 


74.8 ±2.1 




<r 8 (fW0.25) - 47 


0.841 ± 0.076 


0.828 ± 0.070 


0.898 ±0.058 


0.875 ±0.050 


Xeff 


Total 


7533.7 


7540.9 


8065.0 


8074.9 


DIC 


Total 


7557.7 


7565.7 


8089.8 


8097.1 



The models are the same as in Table |nl but here we consider an extended parameter space. We do not fix the primordial helium 
mass fraction, but instead allow it to vary as a free parameter given the condition Y p > 0. We have imposed 7V e ff > 3.046 for 
the cases without sterile neutrinos, and 7V e g > 5.046 for the cases with sterile neutrinos. For the cases without SNe, Ax 2 ff = 7.2 
and ADIC = 8.0. For the cases with SNe from the Union2 compilation, Ax 2 ff = 9.9 and ADIC = 7.3. Further, excluding SN 
data and comparing the sterile neutrino case in an expanded space to the case without sterile neutrinos in a minimal space 
(Table|rrjl, we find Ax 2 ff = 3.0 and ADIC = 10.8. Including Union2 SNe and comparing the sterile neutrino case in an expanded 
space to the case without sterile neutrinos in a minimal space, we find Ax 2 ff = 6.4 and ADIC = 11.9. When adding either N e g 
or f2fc as a single additional parameter to either of the cases with sterile neutrinos, there is a roughly 2a preference above the 
null value. Including Union2 SNe and adding w as a single additional parameter to the sterile neutrino case, there is a 2.7a 
preference for w < —1. This preference is also visible in the analysis with the full extended parameter space shown in this table. 



ated with these datasets. Concretely, the choice of 
light curve fitter in the analysis of SN data dictates 
whether the sterile neutrino model is favored in a com- 
bined analysis of CMB, galaxy power spectrum, Hub- 
ble constant, and SN data. However, the sterile neu- 
trino model is disfavored if we exclude SN measure- 
ments from the analysis. In an attempt to reconcile 
with the laboratory preference for the massive sterile 
neutrino scenario, we proceed to explore if this is an 
indication of new physics beyond the standard cosmo- 
logical model. 

B. Extended Cosmological Parameter Space 
2u s Models 

As discussed in Sec. IIII Al in a combined analysis 
of datasets that include the CMB, galaxy power spec- 
trum, Hubble constant, and SN distances, we obtain 
different results with regard to the viability of sterile 
neutrinos depending on the choice of SN light curve 
fitter. In an analysis without SNe, sterile neutrinos 
are disfavored at A X 2 cff = 11.6 and ADIC = 12.0. We 
examined to what extent this tension could be alle- 
viated in an expanded cosmological parameter space. 
As summarized in Table IIII1 we allow for variations 



in the universal curvature density, constant dark en- 
ergy equation of state, running of the spectral index, 
additional rclativistic species, and primordial helium 
abundance. 

Excluding SN data, adding a single additional pa- 
rameter to the sterile neutrino case does not decrease 
Xoff by a significant amount. For the case of w or 
51^, we find a decrease in xle by about 2, while for 
dn s /dhik, N g, or Y p we find a decrease in xte by 
about 1. For the joint addition of all five of these pa- 
rameters in the model with sterile neutrinos, we find 
that Xefi decreases by 8.6, such that &XeS = 3-0 with 
respect to the ACDM model without sterile neutri- 
nos and no additional parameters. However, due to a 
nonzero Baycsian complexity (see Sec. HI]), we still find 
a large ADIC = 11.5. Hence, including additional pa- 
rameters to the sterile neutrino model decreases ^xls 
to a reasonable level, but the fact that the additional 
parameters are not well constrained is reflected in the 
pessimistic DIC estimates. 

Accounting for the same parameter extension 
(w,£lk, N e g,dn s /dhik,Yp) when adding SN dis- 
tances from the Union2 compilation (i.e. consider- 
ing WMAP+SPT+P(fc)+HST+SNe), we find a de- 
crease in Axgff = 6.4 (down from 17.7) and ADIC = 
11.9 (down from 17.8). When replacing the SNe 
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from Union2 with those from SDSS-SALT2, we find 
Ax^ = 7.4 (down from 20.1) and ADIC = 13.8 
(down from 19.4). Hence, when accounting for SNe 
with the SALT2 fitter, an extended parameter space 
is unable to allow for our 2 massive sterile neutrinos. 

Given the differences between HST and SDSS on 
the best estimate of the Hubble constant @, l68j . 
we considered removing the HST prior on Hq from 
our analysis. We found that excluding the Ho 
prior does not significantly change our constraints, 
mainly because the HST prior only manages to boost 
the best estimate of Ho by about 1 km/s/Mpc 
with respect to the value favored by the CMB and 
large-scale structure data. For instance, consider- 
ing WMAP+SPT+P(fc)+SNe, where the SNe are 
from the Union2 compilation, the Hq constraint lies 
around 70 km/s/Mpc without an HST prior, and 
71 km/s/Mpc when we impose the prior with central 
value around 74 km/s/Mpc. The latter is because 
the data constrains Hq more strongly than the prior 
(such that the error bars on Hq without the prior 
are about 1.4 km/s/Mpc, to be compared with the 
prior of 2.4 km/s/Mpc). This line of reasoning works 
even when excluding SN data. For the particular case 
WMAP+SPT+P(/c), we find A X 2 cS = 9.6 (down from 
11.6) when not including the HST prior. 

We also considered replacing the P(fc) measure- 
ments (with cutoff at k = 0.1 h/Mpc) with two BAO 
distances from SDSS+2dFGRS Considering the 
combination WMAP+SPT+HST+BAO, A x 2 oS = 9.5 
(down from 11.6). Hence, our results are robust to the 
choice of using the power spectrum or BAO distances. 
Moreover, to obtain a better sense of the quoted x 2 g 
values, we note that a universe with w — — 1/3 is dis- 
favored by AXg ff = 96 as compared to a universe with 
w=-l (considering WMAP+SPT+P(fc)+HST). For 
a less extreme case, a universe with w = —0.8 is dis- 
favored by AXcg = 9.4 (with respect to w = — 1). 
These Ax 2 g values significantly increase when further 
including SN data. Hence, our sterile neutrino model 
is disfavored at roughly the same level as a dark en- 
ergy model with w = —0.8 (when not including SN 
data). 

When forcing the 2 sterile neutrinos to be massless, 
Axoff = 5.9 and ADIC = 5.5. Hence, roughly half of 
the degradation in XeS an d DIC could be captured by 
increasing N e s by 2. We note that adding two ster- 
ile neutrinos with a given total mass to = + 7715 
is preferred to adding one sterile neutrino with mass 
to. For example, Ax^ is lower by about 8 when 
{N eS = 5,mi )2 ,3 = 0,to 4 = 0.68 eV,TO 5 = 0.94 eV} 
as compared to {N e g = 4, = 0,m4 = 1.62 eV}. 

However, for a given N e g, the data prefers the sum 
of neutrino masses to be distributed in the least num- 
ber of neutrinos. For instance, given N e g = 5, we 
find that Ax 2 g is lower by about 4 when {7711,2,3,4 = 
0,7715 = 1.62 eV} as compared to {7711,2,3 = 0,TO4 = 
0.68 eV,TO 5 = 0.94 eV}. 

Perhaps more importantly, even when assuming the 



existence of two massive sterile neutrinos, we find a 2a 
preference for an additional massless species. In this 
context, we also note that the extended parameter 
space model with two light sterile neutrinos shows a 
preference for super- acceleration (or w < — 1) [7(j at 
about the 2a level. In fact, this slight preference for 
w < — 1 also persists in a model (with the two light 
sterile neutrinos) that is enlarged only by this one 
parameter (w). 

To summarize, we have studied in detail the ques- 
tion of whether two sterile neutrinos with about eV 
mass each is consistent or disfavored by the latest 
cosmological data. While our sterile neutrino model 
fits each dataset well, in a combined analysis of the 
CMB, Hubble constant, and galaxy power spectrum, 
we have shown that it is difficult to fit all data better 
than a null model without these sterile neutrinos. This 
difficulty persists even when including additional free 
parameters in the cosmological model, such as a con- 
stant dark energy equation of state, curvature of the 
universe, running of the spectral index, effective num- 
ber of neutrinos, and primordial helium abundance. 
Thus, if laboratory experiments continue to favor a 
scenario with two massive sterile neutrinos, and that 
is shown to be at odds with cosmological observations, 
then one may have to look towards a more exotic cos- 
mological model than explored here. 

However, we have also shown that the viability of a 
sterile neutrino model is critically sensitive to our abil- 
ity to identify and control the systematic uncertainties 
associated with the datasets included in our analysis. 
In particular, the sterile neutrino model fits SN data 
better than the null model when using the MLCS light 
curve fitter and worse than the null model when using 
the SALT2 fitter. These differences between the fitters 
can be traced back to different assumptions about the 
nature of color variations in type la SNe and different 
ways in determining model parameters. In a combined 
analysis of CMB, Hubble constant, and galaxy power 
spectrum data, along with SN distance measurements, 
we find that our sterile neutrino model fits the data 
equally well as the null model if we employ the MLCS 
light curve fitter. Thus, a minimally extended model 
with two massive sterile neutrinos could be taken to 
constitute a realistic cosmological scenario, and we ad- 
vocate caution in interpreting combined analyses of 
cosmological datasets given their different systematic 
uncertainties. 



IV. CONCLUSIONS 

Global short-baseline neutrino oscillation data seem 
to favor the existence of two sterile neutrinos with 
masses close to 1 eV each (assuming effectively mass- 
less active species). We have studied the extent to 
which these two neutrinos are allowed by a combina- 
tion of probes including the cosmic microwave back- 
ground, Hubble constant, galaxy power spectrum, and 
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supernova distances. In the analysis of SN data, 
we considered the impact on our results of both the 
SALT2 and MLCS light curve fitters. In particular, 
we showed that the choice of the SN light curve fitting 
method has a major impact on the inferred cosmolog- 
ical model. 

Wc find that the sterile neutrino model provides a 
good fit to each of the considered datasets, and no 
single probe manages to decisively disfavor the ster- 
ile neutrino model with respect to the null model. 
In the joint analysis, sterile neutrinos are allowed by 
the cosmological data (Ax^ ff ps 0) when using the 
MLCS light curve fitter for the SNe in the SDSS 
compilation, and strongly disfavored by the data 
(Axe ff w 18) when using the SALT2 fitter for SNe 
in the Union2 compilation. When excluding the su- 
pernova measurements, the sterile neutrinos are dis- 
favored by the other datasets at Ax^ ff w 12. As an 
illustrative comparison, a cosmological model (with- 
out sterile neutrinos) that has w = —0.8 is disfavored 
by WMAP+SPT+P(/c)+HST (no SN data) at the 
Axpff = 9.4 level compared to the vanilla model with 
w = — 1. 

If the SALT2 fitter is indicative of the correct way 
to interpret SN light curve measurements, then rec- 
onciling two light (~ eV) sterile neutrinos (consistent 
with results from short-baseline neutrino oscillation 
data) with cosmology may require additional freedom 
in the cosmological model. However, no single pa- 
rameter from among nonzero curvature, evolving dark 
energy, additional rclativistic species, running of the 
spectral index, and primordial helium abundance was 
able to decrease AXcfi or ADIC close to zero. In fact, 
even for an extended space with all of these additional 
parameters, the sterile neutrino model is mildly disfa- 



vored at Axe ff ~ 3 (when using the SALT2 fitter). 

The important take-home message, however, is that 
large shifts in Axl s (~ 20) already occur from subtle 
changes to the way parts of the cosmological datasets 
are analyzed. If SN studies converge toward the 
MLCS fitter (as opposed to the SALT2 fitter), then 
two sterile neutrinos with masses close to the eV level 
are easily allowed by the data. Interestingly, even 
when assuming the existence of two massive sterile 
neutrinos, we continue to find about 2a preference for 
an additional massless species. In addition, in this 
model with two sterile neutrinos, a much larger mat- 
ter density would be required (by roughly 40%), which 
helps preserve the constraint on <78(i7 m /0.25) a47 near 
the 0.8-mark, in agreement with galaxy cluster abun- 
dance measurements. The analysis presented in this 
paper shows that it is premature to either rule out the 
existence of two massive sterile neutrinos or claim this 
model is cosmologically preferred. 
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